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We report on the progress we made in studying the infrared behavior of the ghost and gluon 
dressing functions in Landau gauge. Related to this we also investigate a running coupling given 
in terms of those functions and compare our results to those coming from the Dyson-Schwinger 
approach. We present first numerical results for the SU (3) ghost-ghost-gluon vertex renormal- 
ization constant. In addition the spectrum of low-lying eigenvalues and eigenfunctions of the 
Faddeev-Popov operator is determined. The saturation of the ghost propagator in terms of those 
eigenvalues and eigenmodes is discussed at lower momenta. 
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Studying the infrared region in Landau gauge QCD 



1. Introduction 

The predicted infrared behavior of both the gluon and ghost propagator in QCD was the stall- 
ing point of our combined investigation of those Green's functions. In fact, solutions of a truncated 
set of Dyson-Schwinger equations (DSE) for both propagators in Landau gauge predict their dress- 
ing functions Zjj and Zq to follow the power laws ||, f|] 

ZD{q 2 ) « (q 2 ) 2K and Z G {q 2 ) °< ( q 2 )- K with K ss 0.595 (1.1) 

at (extremely) low momentum. Thus, an infrared vanishing gluon propagator occurs in inti- 
mate connection with a diverging ghost propagator. The latter would be in agreement with the 
Zwanziger-Gribov horizon condition as well as with the Kugo-Ojima confinement criterion [|], 0]. 

Furthermore, in terms of the gluon and ghost dressing functions a running coupling a s (q 2 ) 
with a momentum dependence 

a s (q 2 ) = ^Z 2 (q 2 )Z D (q 2 ) with g 2 = 6/j3 for SU(3) (1.2) 

can be defined. This construction is based on the consideration of the ghost-ghost-gluon vertex 
where the renormalization function Z\ of this vertex is assumed to be constant at all momenta. 
Given the asymptotic behavior in Eq. the running coupling would have a finite limit at zero 
momentum 

Those infrared properties are claimed to be independent of how the DSEs are truncated. How- 
ever, it is worthwile to cross-check them using lattice simulations [Qj. We have been pursuing this 
for the SU (3) case during the last two years [|7|]. In addition, we have begun to study the ghost- 
gluon vertex and the spectral properties of the Faddeev-Popov operator. The low-lying eigenmodes 
of the latter are believed to be immediately reflected by a diverging ghost propagator. 

2. The dressing functions and the running coupling 

We have simultaneously studied the gluon and ghost propagators in the quenched approxima- 
tion. Thermalized SU (3) gauge field configurations U = {U x ^ } have been put into the Landau 
gauge. It is well-known that the Landau gauge functional 

1 4 

Mg] = 777LL^Tr^ with g x <=SU(3) (2.1) 

has several local maxima (Gribov copies), each satisfying the lattice Landau gauge condition 
dfiAfi = 0. To explore to what extent this ambiguity has a significant influence on gauge depen- 
dent observables, we have gauge-fixed each thermalized configuration _/V cp times using the over- 
relaxation algorithm, starting from a random gauge copy. For each configuration U, we have 
selected the first (f c) and the best (be) gauge copy (that with the largest functional value) for 
subsequent measurements. For details we refer to our recent publication p]. 

It turns out that the Gribov ambiguity has no systematic influence on the infrared behavior of 
the gluon propagator. In marked contrast to this, the ghost propagator at lower momenta depends 
on the selection of gauge copies. This can be seen in Fig. [I] where in the upper panel the ghost 
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dressing function Zq is presented as a function of momentum q. The lower panel shows the ratio of 
the corresponding functions measured on f c and be gauge copies. We observe that the influence of 
the Gribov ambiguity becomes important at momenta q < 1 GeV. However, the ratio for /3 = 6.0 is 
larger than the corresponding ratio for /3 =5.8 in the sensitive momentum region at the same lattice 
size 24 4 . Therefore, we conclude that the systematic Gribov copy effect at the same momentum is 
reduced if the physical volume is increased. 

We have first studied the infrared behav- 
ior of the gluon and ghost dressing functions. 
In particular, we have tried to confirm the pro- 



3.0 



posed power laws Eq. (1.1). However, fitting 



both functions for q < 0.5 GeV we were un- 
able to extract a common exponent K">0.5.' 
This situation is also reflected in the infrared be- 
havior of the running coupling ct s (q 2 ) which is 
proportional to the product Z 2 } {q 2 )Zo(q 2 ) of the 
two dressing functions (cf. the left hand side of 
Fig. ||). There the fits to the 1-loop and 2-loop 
running coupling are also shown. 

One clearly sees that the running coupling 
monotonously increases with decreasing momen- 
tum as long as q 2 > 0.4 GeV 2 . After passing a 
maximum a s (q 2 ) decreases. One could blame 
finite volume effects for this disagreement with 
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Figure 1: The ghost dressing function Zq mea- 
sured on be gauge copies is shown vs. the (physi- 
cal) momentum q 2 . The lower panel shows the ratio 
(Z fc )/(Z bc ) of dressing functions determined on first 
(f c) and best (be) gauge copies. 



the continuum predictions of the DSE approach. 

Indeed, recent DSE investigations Q on a torus also show a modified infrared behavior of the gluon 
and ghost dressing functions. One could also question whether the ghost-ghost-gluon vertex renor- 
malization function Z\ (q 2 ) is really constant at lower momenta. A recent investigation dedicated to 
the ghost-ghost-gluon- vertex renormalization function Z\{q 2 ) for the case of SU (2) gluodynamics 
[§] supports that Z\{q 2 ) w 1 at least for momenta larger than 1 GeV. We have performed a similar 
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Figure 2: The momentum dependence of the running coupling a s (q 2 ) (left) and of the ghost-ghost-gluon 
vertex renormalization function Z\ (q 2 ) is shown. Both have been measured on be gauge copies. 



'We should mention that the number of data points in this momentum region is certainly too small to fit the asymp- 
totic behavior for q 2 — > 



333/3 



Studying the infrared region in Landau gauge QCD 



study of Z\(q 2 ) for the case of SU(3). Our (preliminary) results are presented on the right hand 
side of Fig. ||. So far the data do not allow to draw a final conclusion whether Z\ (q 2 ) stays constant 
at lower momenta. 



3. Infrared properties of the Faddeev-Popov operator 
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Next we discuss the spectrum of the low- 
lying eigenvalues A; of the Faddeev-Popov (F-P) 
operator. For each maximum of the gauge func- 
tional Fu [g], besides of its discarded eight trivial 
zero modes, the eigenvalues of the F-P operator 
are positive. If all A, > the configuration is 
said to be within the Gribov horizon. 

Having the lowest n eigenvalues and eigen- 
vectors determined, one can construct an ap- 
proximation to the ghost propagator 

G n {k) = T±-® i {k)-® i (-k) (3.1) 

i=\ k i 

with $>i(k) being the Fourier transform of the 
i-th eigenmode. If all n = SV eigenvalues and 
eigenvectors were known, the ghost propagator 
G would be completely determined. In the re- 
cent literature JTo| , 11] it is stated that the F-P 
operator acquires very small eigenvalues in the 
presence of vortex excitations. In any case, ac- 
cording to a popular belief it is an enhanced 
density of eigenvalues near zero which causes 
the ghost propagator to diverge stronger than l/q 2 near zero momentum JjJ]. 

At j3 = 6.2 we have extracted the 200 (50) lowest (non-trival) eigenvalues and their corre- 
sponding eigenfunctions on a 12 4 and 16 4 (24 4 ) lattice. Additionally, 90 eigenvalues have been 
extracted on a 24 4 lattice at j8 = 5.8. This allows us to check that low-lying eigenvalues are shifted 
towards X = as the physical volume is increased. In order to clarify how the choice of gauge 
copy influences the spectrum, the eigenvalues and eigenvectors have been extracted separately on 
our sets of f c and be gauge-fixed configurations. 

In Fig. |3| the distributions of the lowest X\ and second lowest X2 eigenvalues of the F-P operator 
are shown for different (physical) volumes. There h(X,X + A) represents the average number (per 
configuration) of eigenvalues found in the intervall [X,X + A]. Open bars refer to the distribution 
on f c gauge copies while full bars to that on be copies. From this figure it is quite obvious 
that both eigenvalues, X\ and Xi, are shifted to lower values as the physical volume is increased. In 
conjunction the spread of X values is decreased. We have found that the center of those distributions 
tends towards zero stronger than 1/L 2 . Here L refers to the linear lattice extension in physical 



Figure 3: The distribution h{X) of the lowest (left 
panels) and second lowest (right panels) Faddeev- 
Popov eigenvalues Ai and Xi. Filled columns refer 
to be gauge copies, while open columns refer to f c 
copies. 
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units. For example, we have found (Ai)(L) oc l/L 2+e with e sa 0.16(4). It is also visible that the 
eigenvalues X\ and A2 on f c gauge copies are on average lower than those obtained on be copies. 



Eq. ( p.l| ) suggests the popular belief that low-lying eigenvalues and eigenvectors have a major 
impact on the ghost propagator at the lowest momenta. Actually, this is not so easy to predict. We 
have studied, to what extent the lowest n modes saturate the estimator of the ghost propagator G(k) 
on a given set of gauge copies. We show the result in Fig. ^| for the lowest (q\) and the second 
lowest momentum (^2) available on lattices of different size at j3 = 6.2. The estimates for G n (k) 
have been divided by the full (conjugate gradient) estimator G(k) in order to compare the approach 
to saturation for different volumes. 

Considering first the lowest momentum q\ 
we observe from Fig. ^ that the rates of conver- 
gence differ, albeit slightly, for the three differ- 
ent lattice sizes. The relative deficit rises with 
the lattice volume. For the 12 4 and 16 4 lattice 
the rates of convergence are still about the same. 
For example, taking only 20 eigenmodes is def- 
initely too less. One is short of reproducing 
the ghost propagator (by about 50%) whereas 
150 . . . 200 eigenmodes are sufficient to repro- 
duce it within a few percent. For the second 
lowest momentum qi even 200 eigenmodes are 
far from saturating the result. 

In addition, the inverse participation ratio 
has been calculated for the modes in order to 
figure out the localization properties of the F-P 
eigenmodes. It turns out that the lowest eigen- 
modes are not localized on average, although 
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Figure 4: The ratio of the ghost propagator G„(q), 
approximated by the n lowest F-P eigenmodes and 
eigenvalues, to the corresponding full estimator G(q) 
is shown as a function of n for the lowest (q\) and 
second lowest momentum (^2)- The inverse coupling 
is p = 6.2 and the lattice size ranges from 12 4 to 24 4 . 

there are, albeit few, relatively localized ones. 



4. Conclusions 



We have studied the lower momentum region of Landau gauge gluodynamics using Monte 
Carlo simulations on lattice sizes 16 4 , 24 4 and 32 4 . The inverse bare coupling constant has been 
set to /3 = 5.8, 6.0 and 6.2 using the Wilson action. 

Concerning the influence of Gribov copies, it turns out that the gluon propagator does not de- 
pend on the chosen gauge copy while the ghost propagator does. There are indications, however, 
that the strength of this influence, at the same physical momentum, decreases with increasing vol- 
ume. Towards q — > in the infrared momentum region, the gluon dressing function is decreasing, 
while the ghost dressing function is increasing. However, the power laws predicted from the DSE 
approach cannot be confirmed from our data. Correspondingly, the behavior of the running cou- 
pling a s (q 2 ) in a suitable momentum subtraction scheme (based on the ghost-ghost-gluon vertex) 
does not approach the expected finite limit. Instead, this coupling has been found to decrease for 
lower momenta after passing a turnover at q 2 as 0.4 GeV 2 . 
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In addition, we have investigated the spectral properties of the F-P operator and its relation to 
the ghost propagator on a 24 4 lattice at /3 = 5.8 and on 12 4 , 16 4 and 24 4 lattices at j3 = 6.2. As 
expected from another study [Q] we have found that the low-lying eigenvalues are shifted towards 
zero (the configurations approaching the Gribov horizon) as the volume is increased. The low- 
lying eigenvalues extracted on be gauge copies (those with the largest functional value) are larger 
on average than those on f c copies. Thus, for finite volumes better gauge-fixing (in terms of the 
gauge functional) keeps the configurations at larger distance from the Gribov horizon. The study of 
the contributions to the ghost propagator coming from the lowest eigenmodes of the F-P operator 
reveals that even the ghost propagator at the lowest momentum needs not less than 150... 200 
modes for being sufficiently approximated. This result refers to our smallest (12 4 ) lattice at j8 = 6.2. 
For larger volumes and higher lattice momenta the amount of modes needed in an eigenmode 
expansion is even larger. 

All simulations have been done on the IBM pSeries 690 at HLRN. We thank Christian S. Fischer, 
Lorenz von Smekal and Anthony G. Williams for inspiring discussions. We are grateful to Hinnerk Stuben 
for contributing parts of the program code. 
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